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Abstract 

We present a complete algorithm for finding an exact minimal polynomial 
from its approximate value by using an improved parameterized integer re- 
lation construction method. Our result is superior to the existence of error 
controlling on obtaining an exact rational number from its approximation. 
The algorithm is applicable for finding exact minimal polynomial of an al- 
gebraic number by its approximate root. This also enables us to provide 
an efficient method of converting the rational approximation representation 
to the minimal polynomial representation, and devise a simple algorithm to 
factor multivariate polynomials with rational coefficients. 

Compared with the subsistent methods, our method combines advantage 
of high efficiency in numerical computation, and exact, stable results in sym- 
bolic computation, we also discuss some applications to some transcendental 
numbers by approximations. Moreover, the Digits of our algorithm is far less 
than the LLL-lattice basis reduction technique in theory. In this paper, we 
completely implement how to obtain exact results by numerical approximate 
computations. 
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1. Introduction 

Symbolic computations are principally exact and stable. However, they 
have the disadvantage of intermediate expression swell. Numerical approxi- 
mate computations can solve large and complex problems fast, whereas only 
give approximate results. The growing demand for speed, accuracy and relia- 
bility in mathematical computing has accelerated the process of blurring the 
distinction between two areas of research that were previously quite separate. 
Therefore, algorithms that combine ideas from symbolic and numeric com- 
putations have been of increasing interest in the recent two decades. Sym- 
bolic computations are for sake of speed by intermediate use of floating-point 



arithmetic. The work reported in [22l . lid . Il5l . I 111 . |9| studied the recovery of 



approximate value from numerical intermediate results. A somewhat related 
topic is algorithms that obtain the exact factorization of an exact input poly- 
nomial by use of floating point arithmetic in a practically efficient technique 
\^ . In the meantime, symbolic methods are applied in the field of numerical 
computations for ill-conditioned problems 0, H, [23| . The main goal of hybrid 
symbolic-numeric computation is to extend the domain of efficiently solvable 
problems. However, there is a gap between approximate computations and 
exact results jijj]. 

We consider the following question: Suppose we are given an approxi- 
mate root of an unknown polynomial with integral coefficients and a bound 
on the degree and size of the coefficients of the polynomial. Is it possible to 
infer the polynomial and its exact root? The question was raised by Manuel 
Blum in Theoretical Cryptography, and Jingzhong Zhang in Automated Rea- 
soning, respectively. Kannan et al answered the question in [l7j]. However, 
their technique is based on the Lenstra-Lenstra-Lovasz(LLL) lattice reduc- 



tion algorithm, which is quite unstable in numerical computations. In [16 
Just et al presented an algorithm for finding an integer relation on n real 
numbers using the LLL-lattice basis reduction technique, which needed the 
high precision. The function MinimalPolynomial in maple, which finds min- 
imal polynomial for an approximate root, was implemented using the same 
technique. 
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In this paper, we present a new algorithm for finding exact minimal poly- 
nomial and reconstructing the exact root by approximate value. Our algo- 
rithm is based on the improved parameterized integer relation construction 
algorithm PSLQ(r), whose stability admits an efficient implementation with 
lower run times on average than the existing algorithms, and can be used to 
prove that relation bounds obtained from computer runs using it is numer- 
ically accurate. The other function identify in maple , which finds a closed 
form for a decimal approximation of a number, was implemented using the 
integer relation construction algorithm. However, the choice of Digits of ap- 
proximate value is fairly arbitrary [4]. In contrast, we fully analyze numerical 
behavior of an approximate to exact value and give how many Digits of ap- 
proximate value, which can be obtained exact results. The work is regard as 



a further research in 26(]. We solve the problem, which can be described as 
follows: 

Given an approximate value a at arbitrary accuracy of an unknown alge- 
braic number, and we also know the upper bound degree n of the algebraic 
number and an upper bound iV of its height on minimal polynomial in ad- 
vance. The problem will be solved in two steps: First, we discuss how much 
the error e should be, so that we can reconstruct the algebraic number a 
from its approximation a when it holds that \ot — a\ < e. Of course, e is a 
function in n and N. Second, we give an algorithm to compute the minimal 
polynomial of the algebraic number. 

Based on our method, we are able to extend our results with the same 
techniques to transcendental numbers of the form sin -1 (a), log(a), etc., 
where a is algebraic, we also propose a simple polynomial-time algorithm 
to factor multivariate polynomials with rational coefficients, and provide a 
natural, efficient technique to the minimal polynomial representation. The 
basic idea is from [13]. However, we have greatly improved their efficiency. 

The rest of this paper is organized as follows. Section 2 illustrates the 
improved parameterized integer relation construction algorithm. Section 3 
discusses how to reconstruct minimal polynomial and some applications to 
some transcendental numbers by approximations. Section 4 gives some ex- 
perimental results. The final section concludes this paper. 



This paper is the final journal version of 2lj, which contains essentially 
the entire contents of this paper. 
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2. Preliminaries 



In this section, we first give some notations, and a brief introduction on 
integer relation problems. Then an improved parameterized integer relation 
construction algorithm is also reviewed. 

2.1. Notations 

Throughout this paper, Z denotes the set of the integers, Q the set of the 
rationals, R the set of the reals, D(R ra ) the corresponding system of ordinary 
integers, U(n — 1, R) the group of unitary matrices over R, GL(n, O(R)) the 
group of unimodular matrices with entries in the reals, co/jB the i-th column 
of the matrix B. The ring of polynomials with integral coefficients will be 
denoted Z[X]. The content of a polynomial p(X) in Z[X] is the greatest 
common divisor of its coefficients. A polynomial in Z[X] is primitive if its 
content is 1. A polynomial p(X) has degree d if p(X) = ^ d l=0 PiX' 1 with 
Pd 7^ 0. We write deg{p) = d. The length \p\ of p(X) = ^f =0 PiX* is the 
Euclidean length of the vector (po,Pi, • • • ,Pd)] the height \p\oo of p(X) is the 
Loo-norm of the vector(p ,Pi, • • • ,Pd), so Moo = max <j< rf \pi\. An algebraic 
number is a root of a polynomial with integral coefficients. The minimal 
polynomial of an algebraic number a is the irreducible polynomial in Z[X] 
satisfied by a. The minimal polynomial is unique up to units in Z. The 
degree and height of an algebraic number are the degree and height of its 
minimal polynomial, respectively. 

2.2. Integer relation algorithm 

There exists an integer relation amongst the numbers X\,x 2 , ■■• ,x n if 
there are integers ai, a 2 , • • ■ , a n , not all zero, such that Y17=i aiXi = ®- For 
the vector x = [xi, x 2 , ■ ■ ■ , x n ] T , the nonzero vector a = [a 1 , a 2 , ■ ■ ■ , a n ] G Z n 
is an integer relation for x if a • x = 0. 

In order to introduce the integer relation algorithm, we recall some useful 



definitions and theorems 14, 31: 



Definition 2.1. (M x ) Assume x = [xi,x 2 ,--- ,x n ] T G R n has norm |x| = l. 
Define x -1 to be the set of all vectors in R n orthogonal to x. Let 0(R n )nx ± be 
the discrete lattice of integral relations for x. Define M x > to be the 
smallest norm of any relation for x in this lattice. 

Definition 2.2. (H x ) Assume x = [x±,x 2 , ■ ■ ■ ,x n ] T G R n has norm |x| = l. 
Furthermore, suppose that no coordinate entry of x is zero, i.e., Xj 7^ for 
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1 < j < n (otherwise x has an immediate and obvious integral relation). For 
1 < 3 < n define the partial sums 



j<k<n 



Given such a unit vector x, define the n x (n — 1) lower trapezoidal matrix 
#x = (Kj) by 

{0 if 1 < i < j < n - 1, 

Si+i/si if 1 < i = j < n - 1, 

—XiXj/(sjSj + i) if 1 < j < i < n. 

Note that hij is scale invariant. 

Definition 2.3. (Modified Hermite Reduction) Let H be a lower trapezoidal 
matrix, with h^j = if j > % and hjj ^ 0. Set D = I n , define the matrix 
D = (dij) E GL(n, O(R)) recursively as follows: For i from 2 to n, and for j 
from i-1 to l(step -1), set q = nint(hij/hjj); then for k from 1 to j replace 
hi y k by hi^ — qhj^k, and for k from 1 to n replace d^k by — qdj^, where the 
function nint denotes a nearest integer function, e.g., nint(t) = \ t + 1/2J. 

Theorem 2.4. Let 6 R n . Suppose that for any relation m of x and 

for any matrix A G GL(n, O(R)) there exists a unitary matrix Qe U(n-l) 
such that H = AH X Q is lower trapezoidal and all of the diagonal elements of 
H satisfy hjj ^ 0. Then 

1 1 

min — — r < m . 



maxi<j< n _i \hjj\ i<j<n-i \hjj\ 
Proof. See Theorem 1 of bj]. □ 

Remark 2.5. The inequality of Theorem 12.41 offers an increasing lower bound 
on the size of any possible relation. Theorem 12.41 can be used with any 
algorithm that produces GL(n, O(R)) matrices. Any GL(n, O(R)) matrix A 
whatsoever can be put into Theorem 12.41 

Theorem 2.6. Assume real numbers, n > 2, r > 1, 7 > a/4/3 ; and that 
0^x6 R n has O(R) integer relations. Let M x be the least norm of relations 
for x. Then PSLQ(t) will find some integer relation for x in no more than 

'rA /o^(7 n - 1 M x ) 



2 J logr 



iterations. 
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Proof. See Theorem 2 of 14J. □ 

Theorem 2.7. Let M x be the smallest possible norm of any relation for x. 
Let m be any relation found by PSLQ(t). For all 7 > a/4/3 for real vectors 

\m\ < 7 n - 2 M c . 



Proof. See Theorem 3 of [14j. □ 

Remark 2.8. For n=2, Theorem 12 . 71 proves that any relation O^mG 0(R 2 ) 
found has norm |m| = M x . In other words, PSLQ(t) finds a shortest rela- 
tion. For real numbers this corresponds to the case of the Euclidean algo- 
rithm. 

Based on these theorems as above, and if there exists a known error 
controlling e, then an algorithm for obtaining the integer relation can be 
designed as follows: 

Algorithm 2.9. Parameterized Integer Relation Construction 

Input: a vector x, the upper bound N on the height of minimal polynomial, 

and an error e > 0; 

Output: an integer relation m. 

Step 1: Set i := 1, m := 0, r > 2/\/3, and unitize the vector x to 



Step 2: Set by definition 12. 2j 

Step 3: Produce matrix D G GL(n, O(R)) using modified Hermite 
Reduction by definition 12. 3| 

Step 4: Set x := x • ZT 1 , H := D ■ H, A := D ■ A, B := B ■ D~ x , 
case 1: if Xj = 0, then m := coljB; 
case 2: if h^i < e, then m := col n -\B\ 

Step 5: if < |m|oo < N, then goto Step 12; 
if |m|oo > N, there is no such an integer relation, algorithm 
terminating. 

Step 6: i := i + 1; 

Step 7: Choose an integer r, such that T T \h r ^ r \ > T-*\hjj\, for all 
1 < j < n - 1; 
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Step 8: Define a := h rjJ ., j3 := h r+ i }r , A := h r+ i jr+ \, a := a//? 2 + A 2 ; 

Step 9: Change h r to h r+ i, and define the permutation matrix R; 

Step 10: Set x := x • R, H := R ■ H, A := R ■ A, B := B ■ R, if 
i=n-l, then goto Step 4; 

Step 11: Define Q := (qij) G U(n — 1,R), H := H ■ Q, goto Step 
4; 

Step 12: return m. 

By algorithm 12. 9[ we can find the integer relation U(n — 1,R) of the 
vector x = (1, a, a 2 , ■ ■ ■ , a n ). So, we get a nonzero polynomial of degree n, 
which denotes G{x) for the rest of this paper, i.e., 

G(x) = m- (l,x,x 2 ,--- ,x n ) T . (2.1) 

Our main task is to show that polynomial (12. ip is uniquely determined under 
assumptions, and discuss the controlling error e in algorithm 12.91 in the next 
section. 

3. Reconstructing minimal polynomial from its approximation 

In this section, we will solve such a problem: For a given floating number 
a, which is an approximation of unknown algebraic number a, how do we 
obtain the exact value? At first, we state some lemmas as follows: 

Lemma 3.1. Lett be a nonzero polynomial in Z[s] of degree n. If e = 
max!<j< n \a % — di\j, where cti for 1 < i < n are the rational approximations 
to the powers a 1 of an algebraic number a, and a® = 1, then 

|/(«)-/(a)|<5-n-|/U. (3.1) 

Proof. Clear. □ 

Lemma 3.2. Let h and g be two nonzero polynomials in Z[x] of degree n and 
m, respectively, and let a G R be a zero of h with \a\ < 1. If h is irreducible 
and g(a) ^ 0, then 

\g{a)\>n- x -\h\- m -\g\ l - n . (3.2) 



£ is defined by the same way for the rest of this paper. 
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Proof. See Proposition(1.6) of 17|. If \a\ > 1, a simple transform of it 



does. □ 

Corollary 3.3. Let h and g be two nonzero polynomials in Z[x] of degrees 
n and m, respectively, and let a G R be a zero of h with \a\ < 1. If h is 
irreducible and g(a) ^ 0, then 

\g(a)\ > n- 1 ■ (n + 1)"? • (m + 1)^ • \h\^ • \g\^ n . (3.3) 

Proof. First notice that \f\ 2 < (n + 1) • |/|^ holds for any polynomial / of 
degree at most n > 0, so |/| < y/n + 1 • |/|oo- From Lemma [3T2l we get 

\g{a)\ > n- 1 • (n + l)^ ■ (m + 1) V . \ h \^ . ^-n 
This proves Corollary 13.31 □ 

Theorem 3.4. Let a be an approximate value to an unknown algebraic num- 
ber a with degree n > 0, N be the upper bound on the height of minimal 
polynomial of a. For any G(x) in Z[x] with degree n, if 

\G(a)\ < n- 1 ■ (n + l)~ n+1 * • |G|" n • iV 1 "" - n ■ e ■ |GU 

then 

\G(a)\ < n- 1 • (n + 1)^ • |G|" n • N 1 



rl—n 



Proof. Let a G R be with \a\ < 1. From Lemma \3. 11 we notice that \G(a) 
G(a)\ <£• n • \G\oo, and 

\G{a)\ - \G{a)\ < \G{a) - G{a)\, 



so, 



|G(a)| < \G(a)\ +n-e- {G^. (3.4) 
From the assumption of the theorem, since 

\G(a)\ < n- 1 ■ {n + • |G|~" • N^ n -n-e- |GU (3.5) 

combined with (j3.4j) . we have proved Theorem 13. 41 □ 
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Corollary 3.5. Let a be an approximate value to an unknown algebraic 
number a with degree n > 0, N be the upper bound on the height of minimal 
polynomial of a. For any G(x) in Z[x] with degree n, if \G(a)\ < n^ 1 ■ (n + 
l)~ n +h . . N 1 -™, then 

G(a) = 0. (3.6) 

The primitive part of polynomial G{x) is the minimal polynomial of algebraic 
number a. 

Proof. (Proof by Contradiction) Let a G R be with \a\ < 1. According to 
Lemma [3.21 suppose on the contrary that G(a) ^ 0, then 

\G{a)\ > n- 1 ■ (n + l)~ n+ ^ ■ \G\^ n ■ iV 1_r \ 
From the assumption of the corollary, we have 

\G(a)\ < n- 1 ■ {n + l)" n+ i • |G|^ ■ N^ n . 
However, it leads to a contradiction. So, G(a) = 0. 

Let G(x) = Y^ T i=o a i x \ which is constructed by the parameterized integer 
relation construction algorithm from the vector x = (1, a, a 2 , ■ ■ ■ , a n ). Since 
algebraic number a with degree n > 0, according to the definition of minimal 
polynomial, then the primitive polynomial of G(x), denoted by pp(G(x)). 
Hence pp(G(x)) is just irreducible and equal to g(x). Of course, it is unique. 
This proves Corollary 13.51 □ 

3.1. Obtaining minimal polynomial by approximation 

If a is a real number, then by definition a is algebraic exactly, for some 
n, the vector 

(l,a,a 2 ,--- ,a n ) (3.7) 

has an integer relation. The integral coefficients polynomial of lowest degree, 
whose root an algebraic number a is, is determined uniquely up to a con- 
stant multiple; it is called the minimal polynomial for a. Integer relation 
algorithm can be employed to search for minimal polynomial in a straight- 
forward way by simply feeding them the vector (13.71) as their input. Let a 
be an approximate value belonging to an unknown algebraic number a, con- 
sidering the vector v = (1, a, a 2 , ■ ■ ■ , <5 n ), how to obtain the exact minimal 
polynomial from its approximate value? We have the same technique answer 
to the question from the following theorem. 
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Theorem 3.6. Let a be an approximate value belonging to an unknown al- 
gebraic number a of degree n > 0. If 

£ = \ a -a\ < l/(n 2 (n + l) n -5jV 2 ™), (3.8) 

where N is the upper bound on the height of its minimal polynomial, then 
G(a) = 0, and the primitive part of G(x) is its minimal polynomial. 

Proof. Let a G R be with |a| < 1. From Theorem 13.41 and Corollary 13.51 it 
is obvious that 

G(a) = 0, 

if and only if 

\G{a)\ < n- 1 ■ (n + l)~ n+ ^ • |G|~™ • N 1 ' 11 . (3.9) 

Under the assumption of the theorem, we get the upper bound of degree n 
and an approximate value a belonging to an unknown algebraic number a. 
For substituting the approximate value a in G(x), denoted by G(a), there 
are two cases: 

Case 1: G(a) ^ 0, |G(a)| > 0. We have the inequality f)3.5p holds. i.e., 

< n- 1 ■ (n + iy n+1 * ■ \G\^ ■ N x ~ n -n-e- (G^. (3.10) 

Clearly, the inequality (I3.10p satisfies from the condition (13.81) . This proves 
the Case 1. 

Case 2: G(a) = 0. From Lemma I3~TI we have \G(a) — G(a)\ < n ■ e ■ \G\oo, 
hence |G(a)| < n ■ e ■ \G\ao. In order to satisfy condition (13. 9p . we only need 
the following inequality holds, 

n-e- IGU < n- 1 ■ (n + 1)-"+^ • |G|~ n • A^ 1_n . (3.11) 

From theorem 12. 7\ and algorithm 12.91 in Step 5, \G\oo is not more than N. 
Hence we replace \G\oo by N. So the correctness of the inequality (13. lip 
follows from (13. 8p . This proves Theorem 13.61 □ 

It is easiest to appreciate the theorem by seeing how it justifies the fol- 
lowing algorithm for obtaining minimal polynomials from its approximation: 

Algorithm 3.7. Reconstructing Minimal Polynomial 

Input: a floating number (a, n, N) belonging to an unknown algebraic num- 
ber a, i.e., satisfying (13.81) . 
Output: g(x), the minimal polynomial of a. 
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Step 1: Construct the vector v; 



Step 2: Compute e satisfying 



Step 3: Call algorithm 12.91 to find an integer relation w for v\ 
Step 4: Obtain w(x) the corresponding polynomial; 
Step 5: Evaluate the primitive part of w(x) to g(x); 
Step 6: return g(x). 



Theorem 3.8. Algorithm \3. 7| works correctly as specified and uses 0(n(logn+ 
logN) ) binary bit operations, where n and N are the degree and height of its 
minimal polynomial, respectively. 

Proof. Correctness follows from Theorem 13.61 The cost of the algorithm is 
0(n(logn + logN)) binary bit operations obviously. □ 

3.2. Some applications 

In this subsection, we discuss some applications to the practicalities. The 
method of obtaining exact minimal polynomial from an approximate root 
can be extended to the set of complex numbers and many applications in 
computer algebra and science. 

This yields a simple factorization algorithm for multivariate polynomials 
with rational coefficients: We can reduce a multivariate polynomial to a 
bivariate polynomial using the Hilbert irreducibility theorem, the basic idea 
was described in j^] , and then convert a bivariate polynomial to a univariate 



polynomial by substituting a transcendental number in [24| or an algebraic 
number of high degree for one variate in (7|. After this substitution we can get 
an approximate root of the univariate polynomial and use our algorithm to 
find the irreducible polynomial satisfied by the exact root, which must then 
be a factor of the given polynomial. It can find the bivariate polynomial's 
factors, from which the factors of the original multivariate polynomial can 
be recovered using Hensel lifting. This is repeated until all the factors are 
found. 

The other yields an efficient method of converting the rational approx- 
imation representation to the minimal polynomial representation of an al- 
gebraic number. The traditional representation of algebraic numbers is by 
their minimal polynomials [H, 0, [l3, 19]. We now propose an efficient method 
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to the minimal polynomial representation, which only needs an approximate 
value, degree and height of its minimal polynomial, i.e., an ordered triple 
< a,n,N > instead of an algebraic number a, where a is its approximate 
value, and n and N are the degree and height of its minimal polynomial, 
respectively, denoted by < a >=< a, n, N >. It is not hard to see that the 
computations in the representation can be changed to computations in the 
other without loss of efficiency, the rational approximation method is closer 
to the intuitive notion of computation. 

we also discuss some applications to some transcendental numbers by 
using an improved parameterized integer relation construction method. The 
form of these transcendental numbers is sin -1 (a), cos -1 (a), log(a) etc., where 
a is an algebraic number. Moreover, a large number of results were found 
by using integer relation detection algorithm in the course of research on 
multiple sums and quantum field theory in (l2| . 

Suppose (3 is the principle value of sin -1 (a) for some unknown a, which 
is, however, known to the algebraic of degree and height at most n and N, 
respectively. We consider inferring the minimal polynomial of a from an 
approximation to in the deterministic polynomial time. We show that if 
\0 - 0\ is at most s = l/(n 2 (n + l) n "2jV 2n ), this can be done. The specific 
technique is similar with the method in [Tt( ] . 

Thus, in polynomial time we can compute from an approximation a to 
an unknown algebraic number a such that \a — a\ < e, with e as above. Now 
Theorem ( 13.61) guarantees that we can find the minimal polynomial of a in 
polynomial time. 

4. Experimental results 

Our algorithms have been implemented in Maple. The following examples 
run in the platform of Maple 12 under Windows and PIV2.53GB, 512MB of 
main memory. Table [1] proposes the Digits of approximate values to compare 
our method against the LLL-lattice basis reduction algorithm. 

In Table [U we present many examples to compare our new method against 
the LLL-lattice basis reduction algorithm. For each example, we construct 
the irreducible polynomial with random integral coefficients in the range 
—20 < c < 20. Here n and N denote the degree of algebraic number and the 
height of its corresponding minimal polynomial respectively; whereas Dm 
and Dpslq are relative Digits to obtain the exact minimal polynomial in 
theory, E PS lq is in our optimal experimental results respectively. 
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Table 1: Comparison between our algorithm and LLL-lattice basis reduction technique 

From Table [H we observe that the Digits of our algorithm is far less than 
the LLL-lattice basis reduction technique in theory. However, the Digits of 
our algorithm is more than that of the optimal experimental results. So, in 
the further work we would like to consider improving the error controlling. 

The following first two examples illuminate how to obtain an exact quadratic 
algebraic number and minimal polynomial. Example 3 uses a simple exam- 
ple to test our algorithm for factoring primitive polynomials with integral 
coefficients. 

Example 4.1. Let a be an unknown quadratic algebraic number. We only 
know an upper bound of height on its minimal polynomial N = 47. Ac- 
cording to theorem 13. 6[ compute quadratic algebraic number a as follows: 
First obtain control error e = 1/(12 *^*N 4 ) = 1/(1807729447692 * v^) « 
1.0 x 10~ 8 . And then assume that we use some numerical method to get an 
approximation a = 11.937253933, such that \a — a\ < e. Calling algorithm 
13.71 yields as follows: Its minimal polynomial is g(x) = x 2 — 8 * x — 47. So, 
we can obtain the corresponding quadratic algebraic number a = 4 + 3y/7. 

Remark 4.2. The function identify in maple 12 needs Digits:=13, whereas 
our algorithm only needs 9 digits. 

Example 4.3. Let a known floating number a belonging to some algebraic 
number a of degree n = 4, where a = 3.14626436994198, we also know an 
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upper bound of height on its minimal polynomial N = 10. According to 
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theorem 13.61 we can get the error e = l/(n (n + l) n_ 5jV ) = 1/(4 * 5a 
10 s ) ~ 2.2 x 10~ 12 . Calling algorithm 13. 7\ if only the floating number a, such 
that | at — a\ < e, then we can get its minimal polynomial g(x) = x A — 10 * 
x 2 + 1. So, the exact algebraic number a is able to be denoted by < a >=< 
3.14626436994198, 4, 10 >, i.e., < V% + >=< 3.14626436994198, 4, 10 >. 



Remark 4.4. In the example 14.31 we only propose a simple example to repre- 
sent the exact algebraic number by approximate method. In the further work, 
we would like to discuss the efficient arithmetic algorithms for summation, 
multiplication and inverse of the rational approximation representation. 

Example 4.5. This example is an application in factoring primitive polyno- 
mials over integral coefficients. For convenience and space-saving purposes, 
we choose a very simple polynomial as follows: 



V 



3x 9 - 9tr 8 + 3x 7 + 6x 5 - 27x A + 21x 3 + 30x 2 - 21x + 3 



We want to factor the polynomial p via reconstruction of minimal polyno- 
mials over the integers. First, we transform p to a primitive polynomial as 
follows: 
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p = X 



- 3x s + x 7 + 2x 5 - 9x 4 + 7x 3 + 10x 2 - 7x + 1, 



We see the upper bound of coefficients on polynomial p is 10, which has 
relation with an upper bound of coefficients of the factors on the primitive 
polynomial p by Landau-Mignotte bound |20j]. Taking N = 5, n = 2 yields 
e = l/(2 2 * (2 + l) 2 -5 * 5 4 ) = 1/(7500 * \/3) « 8.0 x 10~ 5 . Then we compute 
the approximate root on x. With Maple we get via [fsolveQo = 0, x)]: 
S = [2.618033989, 1.250523220, -.9223475138, .3819660113, .2192284350]. 

According to theorem I3.6[ let a = 2.618033989 be an approximate value 
belonging to some quadratic algebraic number a, calling algorithm 13 . 71 yields 
as follows: 

Pi = x 2 — 3 * x + 1. 
And then we use the polynomial division to get 

p 2 = x 7 + 2 * tr 3 - 3 * x 2 - 4 * x + 1. 



Based on the Eisenstein's Criterion 18], the p 2 is irreducible in Z[A]. So 



the pi and p 2 are the factors of primitive polynomial p. 
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5. Conclusions 

In this paper, we propose a new method for obtaining exact results by 
numerical approximate computations. The key technique of our method is 
based on an improved parameterized integer relation construction algorithm, 
which is able to find an exact relation by the accuracy control. The error e in 
formula (13. 8p is an exponential function in degree and height of its minimal 
polynomial. The result is superior to the existence of error controlling on 
obtaining an exact rational number from its approximation in [26j . Using 
our algorithm, we have succeed in factoring multivariate polynomials with 
rational coefficients and providing an efficient method of converting the ratio- 
nal approximation representation to the minimal polynomial representation. 
Our method can be applied in many aspects, such as proving inequality state- 
ments and equality statements, and computing resultants, etc.. Thus we can 
take fully advantage of approximate methods to solve larger scale symbolic 
computation problems. Furthermore, our basic idea can be generalized easily 
to complex algebraic numbers. 
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